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Abstract. We present a detailed analysis of the morphology and kinematics of nebulae around LBVs and LBV candidates in the 
Large Magellanic Cloud. HST images and high-resolution Echelle Spectra were used to determine the size, shape, brightness, 
and expansion velocities of the LBV nebulae around R 127, R 143, and S 61. For S Dor, R71, R99, and R 84 we discuss the 
possible presence of nebular emission, and derive upper limits for the size and lower limits on the expansion velocities of 
possible nebulae. Including earlier results for the LBV candidates S 119 and Sk— 69°279 we find that in general the nebulae 
around LBVs in the LMC are comparable in size to those found in the Milky Way. The expansion velocities of the LMC 
nebulae, however, are significantly lower — by about a factor of 3 to 4 — than those of galactic nebulae of comparable size. 
Galactic and LMC nebulae show about the same diversity of morphologies, but only in the LMC do we find nebulae with 
outflow. Bipolarity — at least to some degree — is found in nebulae in the LMC as well as in the Milky Way, and manifests a 
much more general feature among LBV nebulae than previously known. 
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1. Introduction 

The most massive stars we know and observe have masses 
above 50 Mq and start as main sequence O stars with lumi- 
nosities of L ~ 10 5 ~ 6 Lq . They inhabit the very upper left part 
of the Hertzsprung-Russell Diagram (HRD). After a hot main 
sequence phase the stars evolve quickly towards cooler temper- 
atures and turn into supergiants within few 10 6 years. Instead of 
evolving further towards the red, the most massive stars enter a 
phase of very high mass loss (about 10~ 4 MQyr -1 ) and reverse 
their evolution back towards hotter effective temperatures. In 
this phase the stars are known as Luminous Blue Variables 
(LBVs). The position of the turning point, and therefore of the 
LBVs (in quiescence) in the HRD, depends on the star's lu- 
minosity, and defines the location of the Humphreys-Davidson 
limit (e.g. Humphreys & Davidson 1979, 1994; Langer et al. 
1994). 

Strong stellar winds and possible giant eruptions in the 
LBV phase peel off more and more of the stellar envelope, and 
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lead to the formation of circumstellar LBV nebulae (e.g. Nota et 
al. 1995; Weis 2001). These LBV nebulae are small, typically 
between 0.2 and about 2pc, and can therefore only be stud- 
ied in our galaxy and, with the higher resolution of the Hubble 
Space Telescope (HST), in a few neighboring galaxies such as 
the Large Magellanic Cloud (LMC). Our knowledge about the 
evolution of the most massive stellar objects in general, and in 
the LBV phase in particular, is sparse. It is not known what 
leads to the very high mass loss rates in these objects and what 
triggers the giant eruptions. These, however, are essential ques- 
tions to understanding the stellar evolution of massive stars. To 
gain insight into the LBV phase and especially the formation 
of LBV nebulae, we studied the nebulae around LBVs in the 
LMC. These LBV nebulae have formed under different con- 
ditions, namely at lower metalicity, than those in the Milky 
Way, and are therefore of great interest. So far the LBV neb- 
ulae in the Magellanic Clouds are the only nebulae, other than 
the galactic ones, which we are able to resolve spatially. In this 
paper we will compare the morphology and kinematics of neb- 
ulae around LBVs in the LMC with those in our Galaxy. 

Only a few LBVs are known (roughly 40, several of which 
are still marked as candidates), of which 10 are in our Galaxy 
and 10 in the LMC. According to Humphreys & Davidson 
(1994), the following objects in the LMC are LBVs: S Dor, 
R 143, R 127, R 1 10, R71, HDE 269582, and LBV candidates: 
S 119, S61, R84, R99. We will add Sk-69°279 to the list 
of candidates (now numbering 11 LBVs/LBV-candidates in 
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the LMC) for reasons discussed in Weis et al. (1997), and 
strengthened in Weis & Duschl (2002), where this object was 
analyzed in detail. Among the LBVs in the LMC, R 143, 
R 127, S 1 19, S 61, and Sk-69°279 are known to possess a cir- 
cumstellar nebula. In this paper we present a study of LBVs 
and LBV candidates in the LMC, concentrating especially on 
the analysis (morphology and kinematics) of the nebulae — 
if present — using high-resolution Echelle spectra and Hubble 
Space Telescope images. The stars R 1 10 and HDE 269582 had 
to be excluded from this work due to a lack of data. Spectra and 
HST images of the LBV candidate S 1 19 have been previously 
analyzed (Weis et al. 2003). The results will be added and put 
into context with this work in the final discussion section. In the 
following we always assume a distance to the LMC of 50 kpc 
(see e.g. Kovacs 2000; Panagia et al. 1991; Westerlund 1990). 

2. Observations and data reduction 

2. 1. Imaging with the HST 

To study the morphology of the quite small nebulae around 
LBVs, images from the HST are especially useful. The fol- 
lowing objects have been observed with the HST (see Tabled 
and analyzed in this paper: R 127, R 143, S61, R71, S Dor, 
R 99. For S 1 19, see Weis et al. (2003). For R 1 10, R 84, HDE 
269582, and Sk-69°279, no HST observations exist so far. 

All objects were observed with the Wide Field Planetary 
Camera 2 (WFPC2). The F656N filter was selected and mim- 
ics quite well an H Q filter, since the radial velocities of the 
LMC stars is roughly 250km s , and therefore H a is within 
the maximum throughput of the filter. All available data were 
retrieved from the STScI data archive and reduced with the 
usual routines in STSDAS/IRAF. All (typically there were 4 
images) longer exposures (500 s) of one object were combined 
and cosmic-ray cleaned. Information about the HST datasets 
are compiled in Tabled The images were not rotated, to ensure 
the full resolution was maintained. The celestial directions are 
therefore indicated in the images. The roll angle of the HST im- 
ages are also given in the last column of Table All HST im- 
ages were taken with the relevant stars centered on the PC chip, 
which has a sampling of 0.0455"/pixel. Most figures shown 
here of the stars and nebulae contain only this PC section. In 
the case of R 143 and R 99, nevertheless, the full mosaiked im- 
ages are shown to discuss the stars and nebulae in context with 
the surrounding ISM. The sampling of those images is lower 
with 0.0996'7pixel. For S Dor, R 71, and R 99 we subtracted 
the Point Spread Function (PSF), which was generated with 
the Tiny Tim (Krist 1995) software for the corresponding posi- 
tions and filters of each star. In all cases we tended to slightly 
oversubtract to make sure that possible residuals are real. 

2.2. Imaging with the 0.9 m-CTIO-telescope 

In addition to retrieving the HST images, we also obtained 
ground-based observations made with the 0.9m-telescope at 
the Cerro Tololo Inter-American Observatory. These images 
were taken with an H Q filter, which contained the [N II] -lines at 
6548 A and 6583 A. The filter was centered on 6563 A, and the 



FWHM was about 75 A. The images were calibrated with bias 
and sky Hatfield frames of the corresponding nights. The see- 
ing was about l'/4, and the nights were photometric. Exposure 
times ranged between 600 s and 900 s for a single image. The 
scale for all images is 0'.'397 per pixel. All images from the 
0.9 m-telescope are displayed with north to the top and east to 
the left. 

2.3. Imaging with the ESO NTT 

We used the ST-ECF/ESO archive to obtain H Q or [N II] images 
of the LBVs, for which no HST images are available, or addi- 
tional information on the larger environment would be useful 
for interpreting our long-slit Echelle data. 

The first archival data set was observed with the red arm 
of the EMMI multimode focal reducer instrument mounted at 
the ESO NTT 3.5 m-telescope. The filter used in these obser- 
vations was an H Q filter (ESO #596) with a central wavelength 
of 6547A and a FWHM of 73 A. The EMMI red arm gives a 
pixel scale of 0'.'27 with the 24 /im pixel CCD (ESO #36). Due 
to limitations from the optical set-up, the usable field size was 
9.'2 x 8'6. Seeing during the observations was 0'.'8, and the ex- 
pose time was 100 s for both fields. We reduced the data in the 
standard manner using IRAF. Cosmic rays were corrected on 
the science frames using the LA-COSMIC IRAF scripts, which 
perform cosmic ray detection using a Laplace filter technique 
(van Dokkum 2001). 

The second archival data set was observed with the SUSI 
imager, also at the ESO NTT. A chronographic unit was in- 
serted, which resulted in a round field with an occulting bar 
running across the field. The set-up and reduction of the data 
are the same as already reported in e.g. Weis (2000). Integration 
time of the R 84 H Q image was 1000 s, and the seeing was 0'.'75. 
The cosmic ray hits were corrected using the LA-COSMIC 
IRAF scripts. 

As for the CTIO 0.9 m data, all ESO NTT images are dis- 
played with north to the top and east to the left. 

2.4. Long -slit Echelle spectroscopy 

High-resolution long-slit Echelle observations were made with 
the Echelle spectrograph on the 4 m-telescope at the Cerro 
Tololo Inter- American Observatory in order to study the kine- 
matics of the nebulae around R 127, R 143 and S 61. Spectra of 
S Dor, R 7 1 , R 84, and R 99 were also taken, in order to search 
for nebular emission possibly connected to the star. All obser- 
vations were made using the same configuration. In this set-up 
for the Echelle spectrograph we replaced the cross-disperser by 
a flat mirror and inserted a post-slit H Q filter (6563/75 A) for 
order selection. The long-slit, which resulted from this config- 
uration, was vignetted to a length of ~ 4'. The spectral region 
that was covered contained the H Q line as well as two [N II] 
lines at 6548 A and 6583 A. We used the 791mm" 1 Echelle 
grating with a slit-width of 150 /im, corresponding to an instru- 
mental FWHM at the H Q line of about 8 km s _1 . The data were 
recorded with the long focus red camera and a 2048 x 2048 
pixel CCD, with a pixel size of 0.08 Apixel -1 along the dis- 
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Table 1. Compilation of the parameters of the datasets 
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HST image 


HST Prog. 


Echelle Spectra 
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R143 


F656N 
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6540 


3 






145.01° 


S Dor 
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S 119 


F656N 


6540 
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148.5° 


Sk-69°279 
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persion and 0'.'26pixel _1 on the spatial axis. The seeing was 
between 1-2" during the observations, and the weather was not 
photometric. Thorium-Argon comparison lamp frames were 
taken for wavelength calibration and geometric distortion cor- 
rection. Telluric lines visible in the spectra were used to im- 
prove the absolute wavelength calibration, the accuracy of 
which is estimated to be 0.04 A (2kms _1 ) or better. 



For each object, the position angle (PA) of the slit was cho- 
sen to match the symmetry of the nebula, as far as the morphol- 
ogy was previously known. Therefore, the PAs are different for 
each object and they are described and shown in images in the 
corresponding sections. The Slit naming was as follows: The 
Slit is named according to the position angle, e.g. Slit PA=222 
is a slit with a position angle of 222 °. If several Slits are taken 
the names indicate the parallel offsets of the slits from the cen- 
ter, e.g., Slit 3N is a Slit 3" north (or S for south) of the Slit 
observed at the stars center. Generally the naming convention 
becomes obvious from figures in each section which show the 
position and names of each Slit for each object. 



The data are presented as echellograms and as measured 
position velocity diagrams (pw-diagrams). All echellograms 
displayed extend 65 A along the spectral axis, which is cen- 
tered on H Q (in rest); in spatial direction the echellograms are 
centered on the projected position of the central star onto the 
slit and are in most cases 1' long. Some echellograms have an 
insert (20" high, 5 A wide) at the upper left, which illustrates 
the emission of the stronger [N II] line (6583 A) again, at differ- 
ent brightness levels to show structures with different surface 
brightness. For R 99 and R 84, the insert includes the total spec- 
trum of the star and not only the [Nil], to better illustrate the 
spectral shape at the star's position in contrast to the emission 
from the larger surroundings. The spectra of S Dor, R99, and 
R 84 are also depicted in their full vignetted length of 4'. The 
pv-diagram of each slit was obtained by measuring the brighter 
[N II] line at 6583 A. The zero position in the pv-diagram cor- 
responds to the projected position of the central star onto the 
slit. All pv-diagrams are generated by binning 4 pixels (1.04") 
along the spatial axis. All velocities are measured in the helio- 
centric system and show a redshift caused by the radial velocity 
of the LMC. 



3. The morphology and kinematics of individual 
objects 

3. 1. R 127 

3.1 .1 . Previous work 

Henize (1955) first recognized that the star R127 (HDE 
269858), in the Large Magellanic Cloud, shows an emission 
line spectrum. Walborn (1977, 1982) classified the star later as 
Ofpe/WN9. The S Dor typ variability (spectrum changed from 
Ofpe/WN9 to an early A, late B) of R 127 led to its classifi- 
cation as an LBV (Stahl et al. 1983). For the stellar param- 
eters of R 127 in the maximum of the S Dor phase, they de- 
rived Mboi = -10.6 m , T cff = 16 000 K, and R, = 150 R©, 
a wind velocity of u w i n d = HOkms -1 and mass loss rate 
of M = 610~ 5 M©yr _1 . The star's mass estimate is about 
60 M© . A good compilation of the photometry and light curve 
of R 127 can be found in van Genderen et al. (1997a). 

Already early low resolution spectra showed typical nebu- 
lar lines, which hinted at the existence of a circumstellar nebula 
(Walborn 1982) and a line split found by Stahl & Wolf (1986b) 
showed the nebula's expansion with roughly 30 — 40kms~ 1 . 
Stahl (1985, 1987) showed, with direct imaging, that the point 
spread function of the star was slightly extended, and added 
evidence for a nebula. Stahl's estimates for the size of the 
nebula was 3" 5 x 4" 5 or 0.8 x 1.1 pc. This measurement was 
supported by long-slit observations (Appenzeller et al. 1987) 
which yield an expansion velocity of v cxp = 28kms _1 , a di- 
ameter of 4", and indicated deviation of the nebula's shape 
from spherical symmetry. Polarimetric observations (Schulte- 
Ladbeck et al. 1993) revealed also that the stellar wind is 
asymmetric. Clampin et al. (1993) were successful in resolving 
the nebula with the use of the John Hopkins Adaptive Optics 
Coronograph. Their H a +[Nll] image with a resolution of Of.'l 
showed a much larger nebula, extending 1.9 x 2.2 pc, and re- 
vealed a diamond shape. In a more recent kinematic study, 
Smith et al. (1998) favored a model with two expanding shells 
around R 127: one inner shell (about 0.6 pc from the star) ex- 
panding with v cxp = 29kms _1 and an outer shell expanding 
with v oxp = 25kms~ 1 . Their abundance analysis made with 
an HST-FOS spectrum at the brighter eastern part of the nebula 
lead to the following parameters of the nebula: log N/H = 8.05, 
log O/H = 8.10, n c ~ 720 cm" 3 , and T c = 6420 K. 
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Fig. 1. This figure shows the HST image of R 127 taken with 
the F656N (H Q ) filter. The field of view is about 15" x 15". A 
north-east vector indicates the celestial orientation. Faint peri- 
odic rings result from the stars point-spread function. 



3.1 .2. The morphology from the HST images 

With the high-resolution HST image, a detailed study of the 
morphology of the nebula around R 127 is possible. Fig. [2 
shows an F656N filter (H a ) image of the nebula. While the cen- 
tral region of the nebula resembles a nearly spherical structure 
with a diameter of 5"4 (1.31 pc), much fainter emission roughly 
north and south adds an elongated shape to the appearance of 
the nebula, henceforth called the Northern and Southern Caps. 
The Northern and Southern Caps extend 1"27 (0.31 pc) and 
l'/72 (0.42 pc), respectively, beyond the central body. The cen- 
tral shell consists of two brighter rims, one at the east, desig- 
nated the East Rim, and one at the (north)-west side, the West 
Rim. The surface brightness is highest at a knot-like structure 
in the West Rim, quite close in projection to a star. The sur- 
face brightness of the nebula is highly variable. Beside obvi- 
ous, distinct low and high surface brightness areas (Caps and 
Rims), knots and filaments of higher surface brightness can be 
found across the entire nebula. These knots reveal sizes of at 
least O'.'l or 0.03 pc for the smallest resolved structures, and 
up to l'/26 (0.3 pc) for the largest (and brightest) knot visible 
in the West Rim. The faintest structures in the nebula are the 
Caps with 2.8 10~ 14 ergscm~ 2 s _1 arcsec -2 for the Southern 
Cap and 3.7 10 -14 ergs cm" 2 s _1 arcscc -2 for the Northern. 
Typical regions in the nebula reach a surface brightness of 
1.2 10 -13 ergscm~ 2 s _1 arcscc -2 , the East Rim is about a fac- 
tor of 2.5 higher, and the high surface brightness knot in the 
West Rim is about 4.5 times brighter than this. 
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Fig. 2. Here, the same image of R 127 as in Figgis shown, 
and the positions of the slits (PA=230°) are marked. Note that 
the slit width was about 1". 



3.1 .3. The Kinematics 

We observed 5 spectra of R 127 with a position angle (PA) 
of PA=230°. The position of the spectra and naming can be 
found in Fig.|2] The corresponding echellograms and position- 
velocity diagrams of each slit for R 127 are shown in Fig. [3] For 
R 127 the top of each echellogram points towards the south- 
west and negative offsets in the pw-diagram are to the north- 
east, positive to the south-west (see also Fig.|2j. 

The radial velocity of the nebula's center of expansion lies 
at about u ra( j = 267km s -1 and is consistent with the stars' 
velocity and H I measurements of the LMC (e.g. Rohlfs et al. 
1984). With an instrumental FWHM of roughly Skms^ 1 , the 
expansion ellipse of the nebula was clearly resolved. At Slit 
Center, which crossed the central star, the expansion velocity 
was 31.5kms -1 . The expansion is much slower in Slit 3N 
and Slit 3S. In Slit 3N the maximum expansion velocity is 
w CX p = 15kms _1 , found at position — 0"5, while the maxi- 
mum expansion found in Slit 3S is v cxp — 12.3 km s -1 at po- 
sition l'.'l. In both slit positions the Doppler ellipse is slightly 
deformed, indicating a deviation from spherical expansion. If 
we determine the spatial center of the ellipse (point of half the 
width) in the pv plots, it would be at — 0"3 in Slit 3N and at 
+ l'/5 in Slit 3S. Both ellipses are much smaller than in the cen- 
tral slit position, and their centers are shifted against each other. 

Futher out, no clear expansion ellipse was detected. 
Nevertheless, the measured velocity values were scat- 
tered around two components — in the case of Slit 6N, at 
243kms~ 1 and260kms- 1 , and for Slit 6S, at267kms _1 and 
283 kms -1 . The emission detected in Slit 6N can be identified 
with the Northern Cap, that in Slit 6S with the Southern. 
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R127 fN JJl 654S H ffi fNlJl65S3 




Fig. 3. Echellograms (left column) and corresponding position-velocity diagrams (right column) of our slits for R 127. 
Measurements are with respect to the heliocentric system. South-west (positive positions) is up, north-east (negative positions) 
down. In the pv diagrams of Slit 3N, Center and 3S models of spherically expanding Doppler ellipses are overplotted. In Slit 
Center several ghost images are present. 
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Fig. 4. This plot shows the stellare spectrum of R 127 as ex- 
tracted from the spectrum in Slit Center. The nebular emis- 
sion in the line of sight to the stars is still superimposed. 
Even though the H a line looks like a P Cygni type profile, the 
blueshifted wing shows additional emission, which is not seen 
in classical P Cygni lines. The spectrum also shows — barely 
recognizable — the three peaks in both [Nil] lines. Next to the 
[N II] line at 6583 A the Fe II at 6587 A is detected. 

Note that at the redshifted side of the Doppler ellipse in Slit 
3N a brighter knot is visible (above the stellar continuum, see 
insert in Fig.[3J, which is identified with the brighter knot found 
in the West Rim. 

In Slit Center an additional velocity component appears 
which has a velocity of 267 km s -1 , similar to the radial ve- 
locity of the star (Stahl & Wolf 1986a) and identical with the 
velocity of the center of expansion. This component is only 
visible at the star's position and does not extend further into 
the nebula. The large spatial size of this component in the pv- 
diagram of Slit Center is only due to the binning of the data 
points. 

The [N II] A6583A/H Q ratio of the nebula around R 127, as 
measured in the spectra reaches 0.6 ± 0.05, this is consistent 
with earlier measurements (see for example Smith et al. 1998). 
The Northern and Southern Caps show the same value. 

A spectrum of the star (within our limited spectral range) 
and the superimposed emission of the nebula is shown in Fig. 
0] The tripel peak in the [N II] emission lines is faintly visible. 

3.1 .4. Discussion 

Combining the morphological information from the HST im- 
ages (Fig. with the kinematics (Doppler shifts) from the 
spectra, it becomes clear that the nebula around R 127 is not 
simply spherical. The nebula consists of an inner, nearly spher- 
ical part, which further out is more and more elongated. This 



elongations is most prominent in the north and south where 
the Caps are attached. The kinematics of the nebula hint that 
even the apparently sherical shell in the center is not expand- 
ing as expected (assuming a constant expansion velocity). The 
deformed expansion ellipses in Slit 3N and Slit 3S indicates 
a deviation from spherical expansion. This becomes obvious 
by comparing the data (asterisks) with the model of an ex- 
pansion ellipse with spherical expansion (line) in the pv plots. 
In both Slit 3N and Slit 3S the data are shifted with respect 
to the sphecial model. If the nebula would be spherical the 
center of expansion should be the same in all slits and lie at 
about 267 km s -1 . Therefore it is concluded that the nebula is 
not expanding with sphercial symmetry in the north and south. 
More precisely, the southern part of the nebula is more red- 
shifed, the northern more blueshifted than the nebula's center. 
The global expansion is bi-directional. This bi-directional ex- 
pansion is supported by the image of the nebula, which shows 
an elogated bipolar shape, most obvious in the Caps. 

In addition to the expanding shell, [N II] emission is also 
present with the radial velocity of the star, but only at the 
star's position. This emission may result either from a sec- 
ond shell which is extremely close to the star, or represents a 
knot-like structure which appears projected onto the star and 
moves more slowly than the shell. The emission might also 
have its origin in a circumstellar disk close to the star which 
was proposed by Schulte-Ladbeck et al. (1993). This third ve- 
locity component (two result from the expansion ellipse of the 
shell) is clearly separated from the first two components of the 
expansion ellipse. It is a single peaked emission line with a 
FWHM=23.6 km s~ ^corrected for the instrumental FWHM). 
A determination from the HST images of a second shell is im- 
possible because of the bright central star. It also shows some 
bleeding which effects the longer exposures and makes it al- 
most impossible to identify any features within 1" of the star. A 
possible second — inner — shell around R 127 was already sug- 
gested by Smith et al. (1998), and they quote an expansion ve- 
locity of 29kms~ 1 . This result is not supported by the new 
measurements presented here and might result from a mixing 
of the three components. 

Most likely the third [N II] emission line structure is in- 
deed a knot within the nebula or closer to the star. This scenario 
would best explain the similarity of the radial velocity of this 
emission line structure with the star's radial velocity. 

3.2. R 143 

3.2.1 . Previous works 

Even though the central cluster of 30 Doradus is home to many 
young and massive stars, only one LBV is known in this region, 
the star R 143 or HD 269929. This star is located roughly 2.'2 
south of R 136, the dense cluster and core of 30 Dor. Based on 
a comparison of different photometric studies with a timespan 
of about 40 years (Parker 1992) and spectral analysis, R 143 
was classified as LBV by Parker et al. (1993). During this time, 
the star changed from a relatively cool star (late F-type) to the 
hot part of the HRD (early B, perhaphs 09.5) and then cooled 
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again to late B. Parker et al. (1993) estimated M bo i ~ -10 m 
and Mzams ~ 60 Mq. 

The star was identified by Feast et al. (1960) and classi- 
fied as an F7 la star. Feast (1961) also identified at least four 
curved filaments close to the star, two of which could belong 
to a nebula around R 143 and are 15" long (3.5 pc, see Figs. 
[5]and|6j. Closer to the star, an elliptical structure was found 
slightly to the west of the central star. The identification of the 
filaments as part of a circumstellar nebula is hindered by many 
filaments and knots which belong to the 30 Dor H II region — 
a clear classification of a circumstellar nebula from morphol- 
ogy alone is not possible. The first high quality images of the 
nebulosity around R 143 were obtained by Clampin et al. (un- 
published, as quoted in Nota et al. 1995). These images show 
stronger [Nil] emission close to the star. Smith et al. (1998) 
analyzed the nebula around R 143 and found that the long fil- 
aments show abundances similar to those of the 30 Dor Hll 
region. Therefore, it was concluded that the filaments are not 
part of an LBV nebula, which should show a higher nitrogen 
content due to the CNO processed material. The filaments also 
move with a nearly constant radial velocity which is of the or- 
der of the velocities found in the 30 Dor complex. However, 
R 143 is surrounded by a much smaller high surface brightness 
LBV nebula. Spectra (Smith et al. 1998) very close to the star 
revealed nitrogen enhanced material. This area coincides with 
the stronger [Nil] emission region seen by Nota et al. (1995, 
Clampins images, see above). About l'/8 north of the star "a 
bona fide LBV nebula" is present. Smith et al. (1998) derive for 
the nebula aT c — 12 200 K, an average n c < 100 cm~ 3 , and 
n e = 1000 cm" 3 in the inner 5". Their [N II] A6583A/H Q ratio 
of the nebula varies according to the distance from the star and 
reaches a maximum of 0.8 at l'/8 north of the star. This neb- 
ula is elongated in the north-south direction with a diameter of 
about 5'.'2 (1.3pc). Their spectra showed a blueshifted motion 
with a velocity difference of 130 km s -1 . They concluded that 
the LBV nebula around R 143 is much smaller, and the fila- 
ments (Feast 1961) are part of the 30 Dor complex and most 
likely not associated with R 143. 

3.2.2. The Morphology from HST images 

The upper left part of Fig. [5] shows an image taken with the 
0.9m-telescope at CTIO (H a +[Nll] filter) and, in the upper 
right, the mosaiked WFPC2 image (F656N), in order to illus- 
trate the position of R 143 with respect to the center of the 
30 Dor complex. The high-resolution WFPC2 images under- 
line the difficulty in disentangling the LBV nebular emission 
from the large background H II region. A section of the PC im- 
age as seen in Fig.[5]includes the LBV nebula. 

The LBV nebula is irregular — triangular shaped — and ori- 
ented north to south-west. The LBV nebula of R 143 shows 
a complex structure consisting of a large number of filaments 
(lower image in Fig.|3 FOV - 10" x 10"). The LBV nebula 
is 4"9 in diameter along the main (long) axis, which corre- 
sponds to a linear size of 1.2pc. The surface brightness is not 
homogeneous, and smaller filaments cross and reach even fur- 
ther out of the nebula, adding a net-like appearance. The rim 
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Fig. 6. A larger section (30" x 30") of the HST image show- 
ing the location of our slits taken from the nebula around R 143, 
including Feasts' filaments. The HST image clearly shows that 
these filaments bend towards the star but are not connected to 
the much smaller LBV nebula. Re-determining the sizes from 
the HST image, the filaments extend to 16"28 (3.9 pc) and 
16"(3.8pc) north-east and north of the star, respectively. 



of the nebula is defined through bent or curved structures. The 
nebula is concentrated to the west of the star, and no coun- 
terpart nebular emission is identified to the east. The surface 
brightness of the LBV nebula close to R 143 was derived to 
~ 110~ 13 ergs cm -2 s _1 arcsec -2 . Feast's filaments are much 
fainter with, on average, 4.3 10 -14 ergscm -2 s -1 arcsec -2 . 



3.2.3. The kinematics 

The nebula around R 143 was observed using three different 
position angles: PA=105°, PA=222° and PA=195°. All slit po- 
sitions are superimposed on the HST image in Fig. [5] 

Using high-resolution long-slit Echelle-spectra, we tried to 
disentangle the emission from the LBV nebula and the back- 
ground 30 Dor H II region using kinematic differences and the 
[N II]A6583A/Hq, ratio — which is expected to be higher in the 
case of an LBV nebula — as indicators. The echelleograms are 
shown in Figs. an d ID as we H as the corresponding pv- 
diagrams. Slit positions PA=195 and Slit PA=222 were used to 
analyze the kinematics of the longer arc-like structures found 
by Feast (1961). In Slit 2N, Center, and 2S the inner tri- 
angular shaped LBV nebula was intercepted. All the echel- 
lograms show a large variety of structures moving with up to 
100 kms -1 velocity differences. This is expected for the re- 
gion around R 143 in the outskirts of the 30 Dor region, which 
shows a complex expansion structure (Chu & Kennicutt 1994). 
In Slit Center the slit crosses the star R 143. Slightly to the 
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R143 




w. 



Fig. 5. This collection of images shows the 30 Dor region (10' x 10') and in subimages the nebula around R 143 first in the 
WFPC2 mosaic and than in the PC of the HST. R 143 is the bright star south-east of the triangular shaped nebula in the small PC 
image. Only this triangular shaped nebula is an LBV nebula. The images illustrated the difficulty to disentangle emission from 
the nebula from that of the 30 Dor H II region. 



north-east of the stars continuum emission (upper direction in 
the echellogram; Fig.0 a [N ll]-bright knot is visible. The pv- 
diagram and the small inset in the echellogram show that this 
knot has two components, one at about 252 km s -1 and one at 
~ 276 km s -1 , a velocity difference of 24 km s -1 . Analogous 
results are found in Slit 2N (where the two velocity components 
are even more obvious) and — though not very prominent — in 
Slit 4N. Slit 2S shows only one component declining in veloc- 
ity from 274 kms -1 to 260 kms -1 . A similar, nearly constant 
velocity component is visible in Slit 4S. In all cases, structures 
with higher velocities are detected, but are not discussed, since 
they can be attributed to the emission from the 30 Dor complex 
and would only lead to confusion if plotted in the pw-diagrams. 



While all the high velocity components in the background 
emission (not included in the pu-diagrams) and the measure- 
ments of the nebula in Slit 4S showed a lower ratio of ~ 0.05 — 
normal for H II regions — the measurements in Slit Center, Slit 
2S and 2N reach values of the [N ll]A6583A/H„ ratio of typ- 
ically 0.5 ± 0.05. The maximum we measured is at 0.7 ± 
0.05. In Slit 4N the [N ll]A6583A/H Q ratio still reaches up to 
0.08. Here due to the seeing, we partly sample the background 
and LBV nebula, which results in the slightly enhanced value. 
The Slits PA=195 and PA=222 cross the star as well as Feasts 
filaments (see Fig. [6ji. The star can be seen in the center in 
both echellograms (Fig. [8} and is surrounded by a small oval 
structure (the actual LBV nebula), which shows at least in Slit 
PA=195 a two velocity component. To the north (downwards 
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Fig. 7. Echellograms (left column) and corresponding position-velocity diagrams (right column) of our slits for R 143 (see also 
Fig. [8}. Measurements are with respect to the heliocentric system. North-west (positive positions) is up, south-east (negative 
positions) down. At Slit Center a ghost image shaped like a dash sign, parallel to the spectral axis, appears, superimposed on the 
[N II]6548 line. 
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Fig. 8. Echellograms (left column) and corresponding position-velocity diagrams (right column) of our slits for R 143. 
Measurements are with respect to the heliocentric system. South-(west) (positive positions) is up, north-(east) (negative posi- 
tions) down. In both images ghost images are on the usual position (see Fig.Qi. 
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Fig. 9. This plot shows a spectrum of the star R 143 with 
a somewhat asymmetric P Cygni profile H Q line and double 
peaked [N II] lines resulting from the LBV nebula in line of 
sight to the star. 



in the echellograms) a longer single velocity structure, slightly 
bended, appears (central emission feature). These emission fea- 
tures (also seen in Slit PA=222) are identified with the fila- 



ments and show no stronger [Nil] emission (see also Smith 
et al. 1998). 

A spectrum extracted from the Echelle spectra of the star 
itself shows (Fig. [9} similarities with the spectrum of R 127. 



3.2.4. Discussion 



Our analysis of the HST images and Echelle spectra support the 
conclusion that Feasts filaments are indeed part of the back- 
ground H II region, and that only a roughly triangular shaped 
high surface brightness nebula close to the star manifests parts 
of an LBV nebula. This LBV nebula shows two velocity com- 
ponents separated by 24kms~ 1 . The small size of the neb- 
ula makes it hard to identify a clear expansion pattern in the 
seeing-limited Echelle spectra. Nevertheless, we conclude that 
the nebula moves with two different velocity components and 
shows a [Nll]A6583A/H Q ratio higher than that of the back- 
ground emission and similar to that of other LBV nebulae. 
The kinematic analysis does not yield other components with 
higher nitrogen emission, which would be part of the LBV neb- 
ula. Smith et al. (1998) find for R 143 a structure that moves 
with a 130 km s -1 difference from the star, which has, accord- 
ing to them, a radial velocity of 348 km s -1 . We find that the 
peak of the stellar H Q line has a radial velocity of 295 kms -1 
and is therefore in better agreement with measurements of the 
Sill at 4128 A, 4131 A, and Mg II at 4481 A by Stahl (pri- 
vate communication), which yield a radial velocity for the star 
of 285.4kms~ 1 , 285.9kms _1 and 287.0kms~ 1 , respectively. 
The large difference of the nebula's and the star's radial veloc- 
ity of 130kms _1 is therefore not supported by our measure- 
ments. 
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The [Nll]A6583A/H a ratio of the LBV nebula we mea- 
sure is comparable to the measurements of Smith et al. (1998), 
with the exception that our maximum ratio is slightly smaller 
(0.7) than theirs (0.8). The nebula's shape and one-sided loca- 
tion (with respect to the star) is very unusual and one can only 
speculate about various scenarios. With R 143 situated in the 
30 Dor region, the strong stellar winds, the larger density, and 
the turbulent motion, which is present in such an H II region, 
exhibit ideal conditions for the easy disruption of circumstellar 
nebulae. 

3.3. S61 

3.3.1 . Previous work 

The star S61 (or Sk— 67°266) was classified as an 08:lafpe 
extr. by Walborn (1977), who noted a similarity of its spec- 
trum with that of R 127. Walborn (1982) puts S 61 in the class 
of Ofpe/WN9 stars and detected double peaked nebular lines 
in the spectrum, indicating that S 61 has an expanding circum- 
stellar nebula (line split 38 km s _1 ). He also noted that nitrogen 
was 13 — 16 times more abundant. Together with the strong 
spectral similarity to R 127 (LBV, see Sect. l3.lt . even though 
no S Dor-type variability (or eruption) is known for S 6 1 , the 
star is a good LBV candidate. Wolf et al. (1987) found that 
the wind velocity of S61 is w W i„d ~ 900kms _1 , and de- 
rived a A/boi = — 10. 3 m . The first observations to resolve the 
nebula around S61 were published by Pasquali et al. (1999). 
The nebula is roughly spherical, being only slightly asymmet- 
ric in the center. North of the star (about l'/2) the nebula is 
brighter, and the diameter is 7'/3 or 1.8 pc. Spectra taken of 
the nebula indicate an expansion velocity, derived from line- 
profile fits, of 28 km s _1 . The nebula and stellar parameters are 
as follows: n c ~ 400 cm" 3 , T c = 6120K (Pasquali et al. 
1999), T* = 36 100 K, R« = 28R , M = 2.2 10~ 5 M yr" 1 
(Pasquali et al. 1997a), while Crowther & Smith (1997) de- 
rive T eS = 27 600 K, i?* = 33 R , M = 1.1 10- 5 M Q yr" 1 , 
r~j 250 kms -1 . 

3.3.2. The Morphology from HST images 

An HST image in the F656N filter is shown in Fig. ITol 
The FOV here is 10" x 10". Previous, ground-based images 
showed that the nebula around S61 is predominantly spheri- 
cal. The new HST image supports this morphology but reveal 
more details of the shape of the nebula. The nebula consists of 
an inner brighter ring-like structure and a larger diffuse emis- 
sion surpassing this ring in all directions. The diameter of the 
inner ring was measured to be 3"4 (0.82 pc) at the smallest 
cross section and 3" 8 (0.89 pc) at the largest . The inner ring 
is not perfectly round, but small dips can be seen. The surface 
brightness is not homogeneous and is brightest in a section of 
the ring to the north. In the ring the surface brightness varies 
around 3 10~ 13 ergscm~ 2 s~ 1 arcsec -2 , while the center of the 
nebula is fainter with about 1 10 _13 ergs cm~ 2 s -1 arcsec -2 . To 
the south and south-west the two largest distortions of the neb- 
ula are visible within the ring. Diffuse emission which shows 
no clear boundary exceeds the inner ring by about 1 " or 0.25 pc. 




Fig. 10. F656N HST image of S 61. The field of view here 
is about 10" x 10". A north-east vector indicates the celestial 
orientation. 
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Fig. 11. Same HST image of S 61 which now shows the loca- 
tion of our Slits. 



Within this more diffuse emission smaller filaments are visi- 
ble. These filaments point roughly in a radial direction away 
from the star. The emission surpassing the ring is slightly 
lower in surface brightness than the nebula's center (about 7 

10~ 14 ergscm~ 2 s _1 arcsec -2 ). 
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3.3.3. The kinematics 

Spectra of S 61 were taken with two different position angles, 
perpendicular to each other — a PA=90° (east-west, 3 spectra 
taken) and PA= 180° (south-north, one spectra taken). The po- 
sitions of the spectra are shown in Fig. ^2 the echellograms 
and position velocity diagrams in Fig.^1 Even though Slit 3N 
is obviously not crossing the nebula, emission is clearly visi- 
ble in this spectrum. Most likely the seeing of 2" is responsible 
for the emission detected at this position. In all echellograms 
(see Fig. I12> either a line split is present, or an asymmetrical 
line shape was found, which indicates a nearly spherical expan- 
sion of the nebula. The maximum expansion velocity was mea- 
sured at position 0" in Slit PA=180 and reaches 26.9kms _1 . 
As expected, the expansion is much smaller in the Slits 3N 
and 2S — 10.5 kms" 1 and 8.9 km s -1 , respectively — which oc- 
curs because the slit intercepts only a section of the sphere. 
Therefore, the diameter of the Doppler ellipse, as well as the 
expansion velocity, is smaller. The shape of the Doppler ellipse 
of our data in Slit Center is slightly asymmetric with a more 
redshifted center of expansion at positive offsets (to the west) 
and a more blueshifted center of expansion at negative offsets. 
This is illustrated by comparing again the data (asterisks) with 
the model of a purely sphercial expansion (lines) in the pv dia- 
grams. Besides this small deviation, the global expansion is in 
good agreement with a spherical expansion. The small devia- 
tion in the expansion ellipse might trace a geometric distortion 
in the line of slight. 

Fig. ^] shows a mixture of an extracted stellar spectrum 
plus the superimposed emission from the nebula which lies in 
the line of sight to the star. The split nebular lines are very 
prominent. The broader wings in the H a emission most likely 
show the contribution of the stellar line. 
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Fig. 13. An extracted stellar spectrum of the star S 6 1 . The split 
lines visible are emission from the superimposed nebula. The 
broad wings seen in H Q manifest the broader component from 
the star itself. 



created at the same time as the ring nebula and expands uni- 
formly and according to the sphercial shell. Finally, the outer, 
faint region of the nebula might be due to an earlier high mass 
loss episode similar to what is found in Planetary Nebulae (e.g. 
Balick & Frank 2002). 



3.3.4. Discussion 

The S 6 1 nebula has the best defined spherical structure among 
the nebulae discussed so far. The HST image shows a ring- 
like morphology, disturbed only by small deformations. The 
same spherical symmetry dominates the expansion pattern of 
the nebula. The nebula around S 61 nevertheless shows differ- 
ences from other classical ring nebulae in being completely sur- 
rounded by fainter diffuse emission, which extends beyond the 
inner ring. The seeing in the ground based spectra and the ex- 
tremely small size of the nebula make it hard to disentangle the 
emission from the ring and the fainter outer part of the neb- 
ula, since both will be observed and sampled in the spectra si- 
multaneously. However, kinematically, the fainter emission can 
not be very much different from that of the ring, otherwise an 
additional fainter component should be visible in the spectra 
with different radial velocity. From the surface brightness de- 
termination of the ring and the fainter outer emission in the 
last section we can estimate that the contribution of the fainter 
emission is about 23%. Whether this fainter, outer emission is 
due to an outflow either in all directions or in the direction that 
is pointed towards the observer (which could form an outflow 
that expands and is seen projected onto the inner nebula) is not 
clear. It might equally be the case that the fainter emission was 



3.4. SDor 

3.4.1 . Previous work 

S Doradus (Sk-69°94, R 88, HDE 35343) is often consid- 
ered the most classical or the prototype of the LBV class. Its 
variability has been known for a long time, and early spec- 
tra showed variability, too (Pickering 1897). Numerous stud- 
ies of its variability have been made, for example by Wesselink 
(1956), Leitherer et al. (1985), Stahl & Wolf (1982), Wolf & 
Stahl (1990), van Genderen et al. (1997b), and Wolf & Kaufer 
(1997). Hubble & Sandage (1953) reported on the brightest 
variable stars in M31 and M33 and already noted, that S Dor 
might belong to this group, later known as the Hubble-Sandage 
variables. Conti (1984) suggested that S Dor type stars — the 
S Dor variables — , the Hubble-Sandage variables, and the so 
called P Cygni variables could be combined into one group, 
now the LBVs. The stellar parameters of S Dor, as measured 
by Leitherer et al. (1985), are T cff = 8 000K, i?» = 300 R , 

M = SlO^MQyr- 1 , and an M bo i 9.7 m taken from 

Lamers et al. (1998). 

Even though S Dor is considered the typical LBV, up to 
now no nebula around the star has been detected. A larger neb- 
ulosity surrounding the star is most likely part of an H II region 
or is a faint superbubble. S Dor has been in its minimum phase 
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Fig. 12. Echellograms (left column) and corresponding position-velocity diagrams (right column) of our slits. Measurements are 
with respect to the heliocentric system. West (positive positions) is up, east (negative positions) down. For Slit PA=180 south 
is up (negativ positions) and north down (positive positions). Slit Center and Slit PA=180 show ghost images. All pv diagrams 
contain the data (asteriks) plus a model of a spherically expanding shell (line). 



during the last few years; but it recently showed an F type spec- 
trum (Massey 1999, 2000). 

3.4.2. Searching for nebular emission 

In Fig. [21 the upper left panel shows an HST image 
(15" x 15") of S Dor in the F656N filter, with (top image) and 
without (bottom panel) the central star. In the bottom left im- 
age the stellar PSF was subtracted using a PSF generated with 



Tiny Tim (Krist 1995), to search for nebular emission close to 
the star. After the PSF subtraction, no nebular emission around 
S Dor was found. The diffuse emission seen between the two 
diffraction spikes in the western direction is most likely not real 
and a result of the scattering light, which produces a ghost im- 
age at that point (see WFPC2 Instrument Handbook, chapter 
5.9), that could not be modeled with Tiny Tim. This residual 
emission is visible in all images of Fig.[H]at the same position. 



14 



K. Weis: Nebulae around LB Vs and LBV candidates in the LMC 



■ 


S Dor 






' M 










SR.. 


.1.. r 


* 



V, 



R71 




• • •••:.;,-/'-vV; 




SDor 



-J 



R71 










R9*9 


■"■ 


■ 










PS 





Fig. 14. HST images of the LBVs S Dor (left panels), R71 (middel panels), and the LBV candidate R99 (right panels). The 
pictures show 15" x 15" sections of the PC field. In the bottom panels the corresponding images are shown, this time with a 
subtracted PSF. 



An image from the 0.9 m-telescope (see Fig. I15t was used 
to study the more distant vicinity of S Dor. A larger, elliptically 
shaped, ring-like nebula surrounds S Dor with a radius along 
the long axis of 5' or about 70 pc. The image shows that this 
structure is clearly defined, and brighter rims are accompanied 
by sections of diffuse emission. To the north-east of the ring 
a brighter H a (+[Nll]) region is visible. The emission of the 
ring around S Dor and the more extended H Q (+[Nll]) emis- 
sion seen in Fig.^l^ 6 psrt of the LMC H II region N 1 19 (or 
DEM L 132). Whether the ring-like structure around S Dor was 
created by S Dor or even resembles an LBV nebula is not clear 
and becomes even more suspicious when comparing the ring 
with a similar ring structure visible at the north-eastern edge of 
the image (Fig. [151. 

3.4.3. The kinematics 

An Echelle spectrum was obtained 3" north of S Dor, running 
east-west. The spatial offset was used to avoid excessive stray 
light and ghost reflections of the grating due to the very bright 
stellar continnum. With the bright central star and a seeing of 
2", stellar emission is still detected in the spectrum (see Fig. 



1 1 61 and could be extracted (see Fig.ll7>. This somewhat noisy 
stellar spectrum, however, shows the star's H Q line with a P 
Cygni profile as well as [N II] lines. The [N II] lines are broad- 
ended with an FWHM (corrected for the instrumental FWHM) 
of about 40km s _1 . The [N II] lines are concentrated on the po- 
sition of the star (see especially Fig. 1161 and are broader than 
in the background [Nil] lines (FWHM ~ 19kms _1 ), which 
are visible in other parts of the spectrum. Since the 40km s _1 - 
broadened emission is visible only at the position of the star, 
they indicate nebular emission very close to S Dor. A determi- 
nation of the [N ll]A6583A/H Q ratio at the star's position is not 
possible since the stellar and nebular H Q emission are super- 
imposed. So far we cannot decide whether or not the detected 
[N II] lines indeed manifest a larger [N II] A6583A/H Q ratio, in- 
dicative of an LBV nebula. 

The echellogram (Fig.1161 is shown in its full length of 4' to 
show the H Q and [N II] emission of the surrounding medium. 
While we found no continuous line split in H Q , which would 
indicate global expansion, the line from the background is also 
broadened (FWHM 28 — 33kms _1 ). This value is compa- 
rable to the turbulent line broadening of an Hll region (e.g 
Chu & Kennicutt 1994). Even though no global expansion pat- 
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Fig. 15. H Q image taken with the CTIO 0.9 m-telescope of the 
larger vicinity around S Dor. The field of view here is about 
10' x 10'. North is up, east to the left. The position of the 
Echelle slit is indicated and drawn to the complete length. 



tern is visible, a second, blueshifted component in the eastern 
(lower) part of the spectrum can be identified. While the broad 
main component of the gas shows emission with a radial ve- 
locity of 281 kms -1 , this second, fainter component (only de- 
tected in H Q ) is moving with 237 km s _1 and is therefore about 
40kms~ 1 slower. This part of the spectrum coincides with the 
brightest section of the H Q emission, visible to the east in Fig. 
IT31 Most likely several layers of filaments in the H II region are 
present here, which move (in part supersonically) with differ- 
ent velocities, as e.g. visible in 30 Dor (see the spectra in the 
field ofR 143 in Fig.|8}. 

Another indication of the origin of the ring around S Dor 
is the [N II] A6583A/H Q ratio. As known from other LBVs, this 
ratio is higher for LBV nebulae than in normal H II regions, 
due to the CNO processed material, assuming that almost no 
mixing with the ISM has taken place since the formation of the 
nebula. The [N n]A6583A/H Q ratio of the ring as measured in 
the Echelle spectra is 0.06, much lower than that of other LBV 
nebulae and comparable to other H II regions in the LMC. If 
the ring was created by S Dor, it was at least not formed during 
the LBV phase. The [Nll]A6583A/H a ratio also makes the in- 
terpretation of the 70 pc diameter ring as faint superbubble un- 
likely, since the typical values are between 0. 1 and 0.2 (Hunter 
1994), due to the diffuse radiation field (and/or low velocity 
shocks) in these objects. 

3.4.4. Discussion 

The HST image of S Dor reveals that, within the limits of these 
observations, no nebular emission is present close to the star. 
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Fig. 16. The echellogram of the spectrum taken 3" north of S 
Dor. The full slit length of 4' is depicted. East is downwards, 
west up. 



The broadened [N II] emission found in our spectrum neverthe- 
less indicate nebular emission very close to and at the star's po- 
sition. Therefore, it might be that the nebula is less than about 
1" or 0.25 pc in diameter. This conservative limit results from 
the PSF subtraction performed on the HST image. Inspecting 
the residuals in the PSF subtracted image, we believe that any 
nebula larger than that should be visible. The limit was mainly 
set by the bleeding of the central star, which prevents a better 
subtraction of the star's PSF. 

On much larger scales, a ring-like nebula surrounds S Dor 
and shows a low [N ll]A6583A/H a ratio. This structure was 
most likely not formed in the star's LBV phase. With an ex- 
pected duration of ~ 25 000 years of the LBV phase (see e.g., 
Humphreys & Davidson 1994), the kinematic parameters — the 
dynamical age is ~ 10 6 yrs — support the arguments that it is 
not an LBV nebula. Nevertheless, we cannot rule out the pos- 
sibility that the ring around S Dor was formed by the star's hot 
wind during the main sequence phase. Weaver et al. (1977) es- 
timated that stars blow a windblown bubble of comparable size 
during their main sequence phase. The ring then is a remnant of 
such a bubble S Dor blew as a hot main sequence star. Whether 
this ring is a main sequence interstellar bubble from S Dor or 
formed purely due to (turbulent) motions within the H II re- 
gion can not be proven from our observations. An LBV nebula 
can nevertheless be ruled out due to kinematics — the dynami- 
cal age is too large — and the chemical composition — no CNO 
processed material is present. 
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Fig. 17. Stellar spectrum of S Dor, extracted from the Echelle 
spectrum. Due to the offset position of the slit, the spectrum is 
of low signal to noise. 



3.5. R71 

3.5.1 . Previous work 

Thackeray (1974) suspected that stars classified Aeq and 
Beq (similar as S Dor in its minimum phase) might indeed 
be S Dor type variables, today's LBVs. The best candidate 
R71 (Thackeray 1974; HD 269006, Sk-71°3) was classified 
B2.5lep (Feast et al. 1960) and has a cool expanding shell 
(Thackeray 1974). Wolf et al. (1981) derived the following 
parameters T c fr = 13 600 K, i?* = 81 Rq, in its minimum 
state, and an Mboi ~ — 8.53 m , rather low for an LBV. For 
the cool expanding envelope they find a temperature of about 
6000K. R71 also shows a strong radiation at 10 /itn, most 
likely due to a cool dust shell (Glass 1984), which is supported 
by IRAS measurements at 12, 25, and 60 /im (Wolf & Zickgraf 
1986). Modeling the dust shell led to the following values: 
Tdust = 140 K, -Rdustshdi = 8000 Rq. In the same year, Stahl 
& Wolf (1986b) detected a broadend, two-component [Nil] 
nebular line in the spectrum of R71, indicating a shell. The 
de-convolution of the line profile into two gaussian, yielded 
a line split of 39kms~ 1 , which would correspond to an ex- 
pansion velocity of about v cxp ~ 20km s _1 . Using ESO 3.6 m 
CASPEC and IUE low resolution spectra, Lennon et al. (1994) 
re-determined the stellar parameters of R 7 1 after re-calibrating 
the extinction curve for LMC mid-B supergiants. Their val- 
ues for R71 are: T cS = 17 250 K, i?„ = 95 R , and a much 
higher, for LBVs more typical A/boi ~ — 9.9 m . Only recently 
was it found that the dust shell of R 7 1 contains amorphous as 
well as crystalline silicate and has a total dust mass of 0.02 Mq 
(Voors et al. 1999). Even though R71 shows the typical light 
variations for an LBV type star (see van Genderen 1979; van 



Genderenet al. 1985, 1988; Lamers et al. 1998) — nebular lines 
and a dust shell — no optical nebula has been detected so far. 

3.5.2. Search for nebular emission 

An HST image of R 71 in the F656N filter is shown in the up- 
per middle panel of Fig. EI as we H as the same HST image 
with the PSF subtracted (bottom panel). After the PSF subtrac- 
tion, no clear indications of an LBV nebula were found. Since 
bleeding affected the images at the star's central position, and 
the PSF model does not include bleeding, residual emission of 
the star is still visible. Testing our PSF subtraction using dif- 
ferently scaled intensities and radii for the PSF models, it was 
concluded that a small, point-like emission to the north-west at 
a distance of 0'.'6 from the star might be real but needs to be 
confirmed. Note that this emission is not the same as the false- 
emission from scattered light reported above, which is more 
arc-shaped and further away from the star. 

3.5.3. A note on the Echelle spectra 

A spectrum (oriented east-west) of the vicinity of R 7 1 was 
taken with an offset of 3" to the north of the star to search for 
nebular emission and expansion. The spatial offset was done 
for the same reasons as discussed for the case of S Dor. The 
spectrum showed no nebular emission at all and is therefore 
not shown here, and does not need to be discussed any further. 
Due to the offset, we definitely missed the possible emission 
knot detected on the HST image. 

3.5.4. Discussions 

From the measurements on the HST images, as well as our 
non-detection of emission in the spectrum, it can be concluded 
that if R 7 1 has an LBV nebula, it is either extremely faint and 
below the detection limit of the HST image and in the spec- 
trum, or very small. If its radius is larger than 0'.'5, it should 
be resolved and visible in the HST images where the PSF was 
subtracted. This limit results again from the reliability of the 
subtracted PSF. At the LMC distance such a nebula would be 
smaller than 0.1 pc. A dust shell around R71 was observed in 
the IR with a radius of roughly 8000 Rq, which would cor- 
respond to 0'.' 00074 (0.0001 8 pc) and would therefore, if emit- 
ting in H Q , not be visible in the HST images. The detection of 
a broadened [Nil] emission (Stahl & Wolf 1986b) is indica- 
tive of an optical nebula, which would expand with less that 
40kms~ 1 . A detection of a point like H Q emission north-west 
of the star (on the PSF subtracted image) might be real and 
could be part of a small nebula around R 7 1 . 

3.6. R99 

3.6.1 . Previous work 

The star R99 (HD 269445, Sk™68°73) was classified as an 
OBf:pe star (Walborn 1977) with P Cygni profiles in the high- 
est Balmer series members. The low-excitation emission spec- 
trum is similar to that of the galactic LBV stars P Cygni and 
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AG Carinae. An IUE study of LMC stars by Hutchings (1980) 
shows that R 99 has unusual high reddening and peculiar emis- 
sion lines. He derived a temperature of T c g = 35 000 K. 
Walborn (1982) and later Bohannan & Walborn (1989) add 
R99 to the list of Ofpe/WN9 stars, a stellar class which in- 
cludes LBVs in their minimum phase (see Stahl et al. 1983, 
1984). 

Analysis of R99 by Walborn (1982) showed (somewhat 
doubtfully) a very weak nebular [Nil] line at 6583 A with a 
line splitting of 61 km s _1 . He notes that this line might be due 
to a circumstellar nebula, but can also be attributed to the sur- 
rounding H II region in which R 99 is embedded. Searching for 
new LBVs, Stahl et al. (1984) combined a large amount of pho- 
tometric data to determine variability and the spectral energy 
distributions of LMC stars. Among them, R99 was found to 
show significant color and brightness variations. They derive an 
extremely high A/boi ~ — 12.2 m . Even though the star shows 
typical LBV characteristics, the authors note that the amplitude 
of the star's variation is small compared to the classical LBVs. 

Stahl & Wolf (1987) note that R99 has a spectral energy 
distribution similar to HD 37836 and speculate that R 99 might 
have a disk. Stahl (1987) searched for a nebula around R99 
using ground-based direct imaging and the subtraction of the 
PSF, without success. Crowther & Smith (1997) again empha- 
sized the peculiarity of the R99's spectrum (showing similar- 
ity with the peculiar LBV candidate HD 5980 in the SMC) and 
its high terminal wind velocity of about lOOOkms -1 . Using 
NLTE atmosphere models and combining NIR, HST, IUE, and 
optical data, Pasquali et al. (1997a, 1997b) obtain the follow- 
ing stellar parameters: T cff ~ 34 000 K, i?„ — 40 Rq, and 
M ~ 3 10~ 5 Mq yr _1 . R99's variability (van Genderen et al. 
1998) supports its LBV membership even though it is quite pe- 
culiar. They suggest that the small variations (see Stahl et al. 
1984) are likely due to a Very Long Term S Dor (VLT SD) 
phase, which shows only a low amplitude. Van Genderen et al. 
(1998) suggest a possible second pulsation mode for R 99. 

Nota et al. (1996a) report the presence of several nebu- 
lar lines, superimposed on the stellar spectrum of R 99, such 
as [Olll], [Nil], and [Sll]. Using the [S II] line ratio, which 
levels at about 1.5, the corresponding density lies at around 
10-100 cm -3 and is therefore close to/at the low density limit. 
Their measurements of the H Q line at different positions across 
the slit identify two major velocity components at ^340 km s _1 
and ^272kms~ 1 . The latter agrees with the centroid of the 
star's H Q line at 271 kms -1 . 

3.6.2. Searching for nebular emission 

The HST image of the star R 99 as in the lower panels in Fig. 1141 
after PSF subtraction does not show any evidence of emission 
from a nearby circumstellar nebula. The image is affected by 
the same artefacts as those of S Dor and R 7 1 ; a ghost image 
shaped like an arc is visible to the north-west. In Fig.^Dan ESO 
NTT-EMMI H Q image is displayed showing a 9' by 8' region 
around R 99. In the same image a blow up shows the closer area 
of R 99 as observed with the HSTs WFPC2 (long edges about 
2.6')- Both images show that R 99 is embedded and surrounded 
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Fig. 19. Echelle spectrum of R 99 and its surroundings, ob- 
served with north-south orientation (south is up). Beside the 
continuum emission of the stars (see upper panel for a closer 
view with different intensity cuts), the spectrum shows H Q and 
[N II] emission of the background H II region. Again, a ghost 
image is present. 



by several H II regions which are identified with the help of the 
plates and catalog of Davies et al. (1976), see caption in Fig. 
n~8l The images illustrate that most of the emission is due to 
the H II regions, and does not show any special centering or 
connection to R99 as it would in the case of a circumstellar 
nebula. 

3.6.3. The kinematics 

We took one spectrum on the star, with the slit oriented north- 
south. The position and length of the spectrum is also shown 
in Fig. CU The echello gram of the full slit is shown in Fig. 
n~9l at the top of which only the stellar spectrum is displayed 
again with different intensity levels. This stellar spectrum was 
extracted and can be seen in Fig. [20] It shows a strong H Q 
emission with two extended and asymmetric wings. The H Q 
line has a very broad FWHM of 157 km s -1 and can be best 
fitted with a Lorentz profile. The H Q emission peak lies at 
286.4 km s~ 1 and therefore differs from the measurents of Nota 
et al. (1996) by 15kms _1 . This difference can be explained 
by taking the lower spectral resolution of their spectra into ac- 
count. Note, also, that taking a rest wavelength of H Q of 6563 A 
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Fig. 18. In this figure we combined images taken with NTT-EMMI and the HST in H Q (+[N n] for NTT-EMMI) of the region 
around R 99. The EMMI image is 9' x 8' large. The footprint of the WFPC2 shows the enlarged region observed with HST. R 99 
is situated in the PC section. The full slit length is indicated in the EMMI image. R 99 lies close to several H II regions in the 
LMC (DEML 158, 169, 166a,b) and is embedded in DEM L 160. Images in the DEM L catalog show even better the large amount 
of H Q emission present at the star's position. Streaks and filaments in the NTT image hint that cooler gas and dust are responsible 
for absorption in these regions. Note that the EMMI image is displayed at a rather soft contrast compared to e.g. the atlas images 
in Davies et al. (1976). 



or 6562.8 A already makes a difference of nearly lOkms -1 . 
Hence, we believe that our measurements are not in contradic- 
tion but an improvement to the radial velocity of the star. A 
component with 340 km s -1 as reported by Nota et al. (1996) 
was not found in our measurements. Clearly, we detect H Q and 
[N II] emission across the entire slit. This is of no surprise since 
we saw that R 99 is embedded in a group of H II regions. The 
peak of the H Q line of the background lies at 294.5 km s -1 
and is therefore 8.1 kms -1 more redshifted than the star. The 
FWHM of the lines, corrected for the instrumental FWHM, are 
28kms~ 1 and 18.5 km s -1 forH Q and [Nil], respectively. The 
[N II] A6583A/H a ratio is about 0.08 ± 0.02 and typical for H n 
regions in the LMC. Walborn (1982) speculated whether the 
faint [N II] emission in his R 99 spectra was from a very faint 
background H II region or a circumstellar nebula. In the same 
manner Nota et al. (1996) discussed that the emission lines 



detected in their spectra can be attributed to a nebula around 
R99. The Echelle spectra show that [N II] is present as part of 
the H II region. The spectrum of R 99 also shows very faintly 
(and somewhat doubtful) the [Nil] at 6584 A (here shifted to 
6590 A). The velocity of this [Nil] line is 294.9 km s" 1 and 
within the errors identical to the radial velocity of the H Q line 
of the background, and rather than the stellar H Q line. We con- 
clude that the [N II] line detected in the stellar spectrum is due 
to the contamination of the spectrum by emission from the 
background and does not represent a nebula around R 99. 

3.6.4. Discussions 

While on a larger scale (see EMMI images in Fig. ED ar, d 
DEML catalog, Davies et al. 1976) R99 is surrounded by 
nebula emission, we find no hints for a circumstellar nebula. 
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Fig. 20. Extracted stellar spectrum of the star R99. The weak 
line at 6590 A is superimposed [N II] background emission (see 
also Fig. fT9i. 

The PSF subtracted HST images in particular lack evidence 
of nebular emission associated with an LBV type nebula of 
R99. Earlier detections of nebular [Nil] emission (Walborn 
1982; Nota et al. 1996) can be accounted by the underlying 
Hll regions — mainly DEML 160 — especially in view of the 
different centroid radial velocities of the nebular lines with the 
stellar H Q line. We therefore conclude that the nebular emission 
detected in spectra of R 99 results from background emission 
which was acquired at the same time of these observations. The 
densities derived by Nota et al. (1996) of 10 — 100 cm ~ 3 are 
more typical values for H II regions and support the suggestion, 
that the nebular emission is part of the background H II rather 
than a circumstellar LBV nebula. If that is the case, the den- 
sity would generally have higher values (several 10 2 ~ 4 cm -3 ). 
The star, however, has an unusally shaped H Q line (asymmet- 
ric Lorentz profile) which is quite broad, again supporting its 
exceptional state in our sample. 

3.7. R84 

3.7 A. Previous work 

R 84, also known as HD 269227 or Sk-69° 79, was classified 
as spectral type Olafpe by Walborn (1977), as he noticed "sharp 
high-excitation emission", P Cygni profiles, and strong hy- 
drogen emission. Cowley & Hutchings (1978) and Hutchings 
(1980) noted that the spectrum of R 84 shows indications for a 
cool supergiant and concluded that the spectrum is composite 
with a BOe and an M2 supergiant, showing clear TiO-bands. 
The conclusion that the spectrum of R 84 is a combination of 
a luminous late type and a hot early type star led to the idea 
that R 84 has a companion star. Later R 84 was included in the 



Ofpe/WN9 sample (see Walborn 1982; Bohannan & Walborn 
1989). R84 showed strong [Nil] nebular lines which accord- 
ing to Walborn (1982) are "quite strong. ..and although there is 
no clear velocity structure, they are broader than single, unre- 
solved nebular lines would be". 

Glass (1984) reports a strong IR excess but finds no TiO- 
bands; he identifies R 84 as a composite Wolf-Rayet star plus 
late type supergiant. Studies by Stahl et al. (1984) confirm the 
large IR excess, but note that the (K-L) value is too high to 
result only from a late-type companion. They propose that a 
circumstellar dust shell surrounds R 84 and strengthen this ar- 
gument with the peculiar variation of the shape of the contin- 
uum of the star and a variable 2200 A-feature, which could be 
due to dust formation. Stahl et al. (1984) derived the following 
stellar parameters for R84: T cff = 25 000 K, R« = 33 R Q , 
and Mboi ~ — 9.3 m . R84 is accordingly described as an S Dor 
variable with small amplitude variations. 

Further evidence for a circumstellar shell around R 84 came 
from a broad, most likely non-interstellar Nal D line and the re- 
solved (FWHM= 37km s" 1 ) [Nil] lines (Stahl & Wolf 1986b). 
Direct imaging observations to resolve a circumstellar nebula 
were not successful (see Stahl 1987). Wolf et al. (1987) raised 
doubts about the fact that the companion star is responsible for 
the star's peculiarities, showing that the UV spectra of R 84 is 
that of an 09.5lab star and quite similar to that of S61 (see 
Sect. I3.3I >. Several atmospheric analyses of R 84 (Schmutz et 
al. 1991; Crowther et al. 1995; Nota et al. 1996a; Pasquali et 
al. 1997a) in the few last years derived consistent results for the 
stellar parameters, also comparable to those mentioned above, 
but favored different evolutionary scenarios for R 84. Schmutz 
et al. (1991) proposed that R 84 had gone through a red super- 
giant phase, Crowther et al. (1995) compared it with the dor- 
mant LBV R71, and Pasquali et al. (1997a) identified R84 as 
a quiescent LBV. 

A recent analysis with ESO NTT SUSI images and using 
ADONIS adaptive optics by Heydari-Malayeri et al. (1997) 
showed that, if a companion star is in the line of sight to R 84 
(or a close binary), it must be closer than (X'12. Van Genderen 
& Sterken (1999) concluded that R84 shows no signs of vari- 
ability and earlier reports of variations might be due to faint 
field stars in the photometric aperture. They classify R 84 as an 
ex-/dormant LBV. Obviously, the question of R 84's evolution- 
ary state and history is not settled yet. We include R 84 in our 
set of LBV candidates, and discuss the possibility that R 84 is 
an LBV and whether it is surrounded by a circumstellar shell. 
No high-resolution or HST images of R 84 are available at this 
point. 

3.7.2. Images 

Since no HST images are available, we used ESO NTT-SUSI 
archive images with a seeing of Of! 7 5 and inspected these im- 
ages for possible nebular emission associated with the star. A 
NTT-SUSI image of R 84 is shown in Fig.|^ together with an 
image on a larger scale (9 ' x 8') taken with the NTT-EMMI. 
In both images an H a (plus [N II]) filter was used. The EMMI 
image (Fig. 12 lb shows that R 84 is in the vicinity of the H II 
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Fig. 21. This figure shows a composite H a image constructed from ESO NTT-EMMI and higher resolution ESO NTT-SUSI 
observations of the area around R84. The EMMI image (left) is 9' x 8', the inlet with the SUSI image is somewhat larger than 
50". The SUSI image was taken using a coronographic mask, which blocks part of the emission in the field. R84 is the bright 
star in the northern part of the SUSI image. Due to its brightness, charge transfer effects are seen starting at the star's position 
and extending down the chip as a long streak. In addition, the slit position with the full 4' length is indicated in the EMMI image. 
The larger H II region north-east of R 84 is part of DEM L 1 13. Faint filamentary nebular emission (DEM L 1 10) is also seen west 
of R 84, stretching roughly north-south. The high density of stars visible particularly in the EMMI image is due to the location 
of the field near to the stellar bar of the LMC. 



region DEM L 1 13 (Davies et al. 1976), which is partially seen 
in the upper left corner of the EMMI image. To the west (right) 
of R 84 very faint H Q filaments are visible, which are part of 
DEML 110 and are roughly oriented north-south before they 
kink to the east. Comparing the EMMI images with the im- 
ages of the DEM L catalog shows that very faint emission also 
lies between the the filaments (DEML 110) and DEML 113 
and would hence be present at the position of R 84. In the ESO 
SUSI image, which shows an area of about 13 x 13 pc, no obvi- 
ous H Q emission was found. R 84, which is the brightest star in 
the northern part of the image, nevertheless show no perfectly 
round structure. Some streaks point away from the star, most 
of them due to diffraction spikes and charge transfer problems 
(long line streching down). 



3.7.3. The kinematics 

We took one Echelle spectrum centered on the star (see Fig.l22i 
and one offset by 3" north, both with an east-west orientation. 
Since the offset spectrum does not contain any additional in- 
formation, it is not shown and will not be discussed here. Of 
the LBV candidates in the LMC, R 84 is the strangest object. 
In its spectrum the H a line shows a broad and a narrow compo- 
nent (see extracted spectrum in Fig.l23i. The broad component 
has an FWHM of 435 km s -1 . The narrow component is split, 
as are the [N II] lines. The peaks of the split have radial ve- 
locities of 241 kms -1 and 265 km s -1 , and are therefore sepa- 
rated by 24kms~ 1 (in agreement with measurements by Stahl 
& Wolf 1986a). The narrow split H Q line seems superimposed 
on the broad component. The narrow and broad component can 
be fitted (and subtracted) with just two gaussians. The split 
[N II] emission is quite bright, and if we compare it with the 
narrow H Q line, we obtain an [Nll]A6583A/H a ratio of about 
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Fig. 22. The Echelle spectrum of the star R 84 and its vicinity, 
slit orientation is east-west. West is up. The upper panel illus- 
trates the stellar continuum spectra in which the broad H a and 
[N II] lines show up at higher intensity. Across the slit some 
background H Q emission was found, but the background [N II] 
is barely visible. 



0.5 ± 0.1. It is probable that different origins are responsible 
for the broad and the narrow components. The narrow compo- 
nent most likely results from nebular emission, while the broad 
component represents H Q emission from the star. The echel- 
logram (see Fig. I22> shows that the nitrogen emission is de- 
tected only at the position of the star (compare full spectrum 
with insert to the top, which contains only the star, but with 
different brightness cuts). While H Q emission of a faint back- 
ground Hll region is seen across the entire slit, [Nil] emis- 
sion is not, or is only barely visible. It is at least much fainter 
than the [N II] emission seen superimposed on the star's con- 
tinuum. Summing up all the [N II] emission for a better signal 
to noise, we obtain an [N ll]A6583A/H Q ratio of 0.1 for the 
background, much less than the ratio at the star's position. In 
the echellogram we can see that the H Q line is slightly brighter 
closer to the star (up and down in the echellogram, Fig. I22> . 
This brighter region has a diameter of 23". At the same time, 
no bright [Nil] is seen, and the [N ll]A6583A/H a ratio stays 
the same. 
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Fig. 23. An extracted spectrum of the star R 84. The nebula 
[N II] lines are doubled as is the narrow H a line. 

3.7.4. Discussions 

We find that bright [N II] emission is concentrated at the po- 
sition of the star and supports the presence of a very small 
circumstellar nebula around R84 which is not spatially re- 
solved. The size of such a nebula can only be estimated to be 
smaller than 2" (about 0.3 pc). This limit is given by the see- 
ing we had in the Echelle observations. The NTT-SUSI image 
shows some faint protrusions which could be part of a small 
circumstellar nebula, but the quality of the data precludes a def- 
inite statement on details of the nebula structure. The spectrum 
indicates that the expansion velocity of this nebula is about 
12kms _1 , as derived from the separation of the narrow lines. 
The [N ll]A6583A/H Q ratio is in good agreement with values 
found for LBV nebulae. Faint nebular emission was also ob- 
served above or below the star in the spectrum but shows a 
much lower [N ll]A6583A/H Q ratio, typical for an H II region. 
The broad component in the spectrum of R 84 — only visible 
in H Q — has not been seen in other LBVs and opens questions 
about whether it is a new spectral feature unique to R84 or 
could possibly be related to a binary. The narrow H Q compo- 
nent is most likely from the nebula that also forms the [N II] 
lines. The broad H Q component could be the stellar emission 
broadened due to a very strong stellar wind. 

4. Summary of previous results for other LBV 
candidates in the LMC 

Observations and analysis for the two LBV candidates S 119 
and Sk— 69°279 have been published earlier in Weis et al. 
(2003) and Weis & Duschl (2002), respectively. In the follow- 
ing we will briefly summarize these results to put these objects 
in context with the results of this paper. 
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Note in this context that we use the term outflow from a 
nebula in the sense of outflowing material from a disrupted 
shell. Outflow here is used to describe the gas which flows out 
of the nebula, it is not used for the stellar wind outflowing from 
the star. 

4.1. S119 

S 1 19 (Sk-69° 17, HDE 269687) is classified as an Ofpe/WN9 
star (Bohannan & Walborn 1989), and the latest determina- 
tion of the stellar parameters yield T cS = 26 200/27 000 K, 
logL = 5.76/5.8OL ,andlogAf = -4.87/ - 4.92M yr- 1 
(Crowther & Smith 1997). Different values are obtained for 
two different models to take the nebula contamination into ac- 
count. The terminal wind velocity is 400km s _1 (Pasquali et al. 
1997a). Nota et al. (1994) were the first to resolve a 7" x 9" 
large nebula (1.9pc x 2.1pc) associated with S 119. We anal- 
ysed this nebula (Weis et al. 2003) using HST imaging and 4 m 
long-slit high-resolution spectra, in the same set-up as descibed 
in this work. The images show that the nebula is nearly spher- 
ical with a diameter of 7'.'5 corresponding to 1.8 pc. The shell 
is brightest in the east. From the south-east toward the north- 
west part of the nebula several filaments are stretching out, the 
longest extending 1'.'81 beyond the nebula's shell. Our long 
slit observations show an expansion of the nebula main shell 
with a maximum expansion velocity of 25.5 km s -1 . Beside 
the spherical expansion of the central part of the nebula, we 
showed that — as the filaments in the images indicate — material 
is streaming out with much higher radial velocity. The highest 
velocity detected in the outstreaming gas is with 283kms _1 
about 130kms _1 faster than the center of expansion (with 
156 km s -1 ). The nebula around S 1 19 is not completely closed 
and shows outflow. Such an outflow might results, e.g., from 
Rayleigh-Taylor instabilities in the nebula, density gradients in 
the ambient medium or the onset of an asymmetric faster stel- 
lar wind. Finally, note that the location of S 1 19 in the disk of 
the LMC is still under debate since the star's radial velocity 
is much slower than expected (Nota et al. 1994; Danforth & 
Chu 2001; Weis et al. 2003). Hence, the nebula's size might be 
different, giving a possibly shorter distance. 

4.2. S/c-69°279 

Sk— 69°279 was recognized as an H Q emission-line star by 
Bohannan & Epps (1974) and later classified as 09f by Conti 
et al. (1986). Thompson et al. (1982) derived an effective tem- 
perature T c g = 30 300 K, an absolute bolometric magnitude 
M bol = -9™72, and a visual one M v = -6™77. We found 
(Weis et al. 1995) that the star Sk-69°279 is surrounded by a 
nebula with a diameter of 18" or 4.5 pc. The nebula expands 
with roughly 14kms _1 . Several analysis and studies of this 
object (Weis et al. 1995, 1997b; Weis & Duschl 2002) let us 
suspect that this object is an LBV candidate due to its bright- 
ness, UV spectra (see Smith Neubig & Bruhweiler 1999), and 
the large nebula, which is especially strong in the [N II] lines. 
The [Nll]A6583A/H Q ratio of 0.7 is comparable to those of 
other LBV nebulae (e.g., AG Car: [N ll]A6583A/H Q - 0.5, de 



Freitas Pacheco et al. 1992). Recently (Weis & Duschl 2002), 
we detected that Sk— 69°279 also exhibits a large scale out- 
flow similar to that seen in S 119. This outflow is seen as a 
7'.'1 (1.7pc) large filament extending to the north of the nebula. 
The filament moves about 21 kms -1 faster than the center of 
expansion of the nebula. 

5. Summary of this work 

From all other stars analyzed here, R99 was the object that 
showed the least evidence of a nebula. Neither the PSF sub- 
tracted HST images, nor the Echelle spectra show any indica- 
tions of nebular emission. If an LBV nebula surrounds R 99, 
it must be at least smaller or much fainter than what we can 
measure, giving our detection limits. The LBV/LBV candidates 
R 127, R 143, S61, S 119, and Sk-69°279 all show nebulae 
which are several arcseconds in diameter. We find that R 127's 
nebula is fairly spherical with two triangular-like attachments, 
the caps. The kinematics of the nebula gives hint for bipolar- 
ity. In R 143 only a smaller irregularly — nevertheless approxi- 
mately triangularly — shaped nebula with net-like structures is 
identified as an LBV nebula. Our spectra show two velocity 
components, separated by 24kms _1 associated with this LBV 
nebula. For S61, the nebula's morphology is predominantly 
spherical (~3'.'6), with fainter emission that — in projection — 
surpasses the borders, the global expansion velocity is about 
27kms~ 1 . Analyzing images of S Dor, R71, and R84, we 
found no clear indication for nebular emission. Filamentary, in 
the case of S Dor even bubble-like emission visible in H Q found 
close to S Dor and R 84 are part of the H II regions or diffuse 
hot gas in the LMC, but are not relics of the stars' LBV phase. 
However, note that all these objects, superimposed on the stel- 
lar spectrum, show nebula [N II] lines which are broadened (S 
Dor) or split (R71, R 84). These lines are indicative of nebular 
emission in connection with these stars. The higher line ratio 
and spatial extend might result from LBV nebulae close to the 
stars which we can not resolve spatially. 

6. General Discussion and Conclusions 

In this final section we will compare the parameters of the LMC 
LBV nebulae studied here with those of the known Galactic 
LBVs. Therefore, we summarize the morphologic and kine- 
matic parameters of both groups in Table [2] For a detailed de- 
scription and sources of the parameters, we give examples of 
the most recent literature which, however, is far from complete. 
In Table|2]we quote the size, expansion velocity, and morphol- 
ogy of the nebulae. In some cases, where no clear expansion 
ellipse could be traced, not the expansion velocity but the line 
split is given, which is marked with (split). For S Dor the line 
was not split, so the FWHM is given instead. 

6. 1. Morphology and Nebula Sizes 

The different morphologies and sizes of nebulae around LBVs 
and LBV candidates in the LMC are particularly conspicuous. 
The nebulae are spherical, bipolar, irregular, or show outflow. 
Also, they appear to be grouped into two different size classes. 
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Either they are quite large ( 1 pc or larger) or they are very small 
(upper limit from the resolution of our data 0.3 pc). Of course, 
we have to keep in mind that selection effects play an important 
role. In the LMC we predominantly detect the largest nebulae 
since those smaller than about 0.5 pc correspond to an appar- 
ent size of about 2" in the LMC and are therefore much harder 
to find in surveys, for instance. This is especially the case for 
stars which are not identified as LBVs so far because they are 
in a dormant state. The LBV nebulae in the LMC are in gen- 
eral about the same sizes (perhaps slightly larger) as those in 
our Galaxy. From all nebulae which have been resolved so far 
this shows that they cover sizes between 0.15 (HD 168625) 
and 4.5 pc (Sk— 69°279). Taking again all resolved nebulae into 
account, we can estimate an average size for LBV nebulae of 
~1.25pc. 

LBV nebulae in the Galaxy and the LMC come in the same 
variations of morphologies. In particular bipolarity — at least to 
some degree — in LBV nebulae is present in nebulae in both 
galaxies. Interestingly, we find that bipolarity is found in neb- 
ulae of very different sizes. A bipolar nebula par excellence is 
seen in the Homunculus around rj Car, but also larger nebulae 
like those around HR Car and R 127 in the LMC have bipo- 
lar components. Conclusively, bipolarity is neither restricted to 
size, nor connected to the host galaxy (or as a consequence: 
metalicity). 



6.2. Expansion velocities 

While the sizes are about the same for LBV nebulae, the ex- 
pansion velocities are significantly different. Looking at nebu- 
lae in our Galaxy, it seems at first glance the larger a nebula, 
the lower its expansion (Fig. [24}. This can manifest some kind 
of evolution — in connection to a slow-down process — of LBV 
nebulae. As argued in Weis et al. (1997a) and Weis (2001), HR 
Car does show signs of being an aged, slowed down version of 
r\ Car's nebula. We would expect that larger nebulae are older 
since they had more time to expand, but at the same time the 
expansion will have slowed them down. However, this scenario 
of larger nebulae being slower completely breaks down for the 
nebulae in the LMC. First, we note that all LBV nebulae in the 
LMC are expanding much more slowly than those in the Milky 
Way. In the LMC none is moving faster than 32kms~ 1 — in the 
Milky Way the expansion velocities are up to an order of mag- 
nitudes higher. The slowest LBV nebula in the Milky Way is 
about as fast as the fastest LMC nebula. As a consequence, the 
spread in expansion velocity of nebulae in the LMC is much 
smaller than in the Galaxy. While we find velocities between 
26kms~ 1 and several 1000 km s -1 in the Galaxy, the LMC 
LBVs range from only 12 to 32kms~ 1 . The expansion veloci- 
ties of LBV nebulae in the LMC are confined to a small range. 
From the presently known sample we conclude that there is a 
significant difference in the expansion velocities between the 
nebulae in the LMC and the Milky Way. Even if we exclude 
r\ Car at that point, given its recent outburst and exceptional 
status, this statement holds true. Note in this context, that de- 
tecting higher expansion velocities in nebulae around LBVs in 
the LMC is not a problem. Even if we are not resolving the 
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Fig. 24. In this plot we compare the expansion velocities and 
sizes of nebulae around LBVs in the Milky Way (triangles) and 
the LMC (circles). A dashed line indicates the limit (0.3 pc) up 
to which we can resolve nebulae in the LMC. In the case of 
asymetric nebular lines show the limits in each symetry axis or 
connect datapoints of the same object in case two nebula parts 
exists (P Cyg and -q Car). 

nebulae, the high expansion velocity will be visible in spectra 
which include the [N II], which will be broadend, splitted or at 
least be shifted against the stellar radial velocity. 

So what determines the expansion velocities of LBV neb- 
ulae, and why are they so different between the LMC and the 
Milky Way? A lot of different stellar and environmental pro- 
cesses possibly affect the expansion velocities, only the most 
obvious of which we will discuss in the following. Most likely, 
however, a mixture of the various processes will occur. 

- Stellar winds: The mass loss history and change of wind 
velocity during the LBV phase is a function of the star's 
initial mass, for instance. The formation of an LBV neb- 
ula is thought to occur due to wind-wind interaction, as 
proposed in, e.g., Garcfa-Segura et al. (1996). If, now, the 
input for wind velocity and mass loss is different, the ex- 
pansion velocity of the nebulae will naturally differ due to 
the solution of the wind-wind interaction models. The ex- 
pansion velocities will also change, if the duration of the 
LBV phase (or more precisely the phase of slow wind) is 
different, since then the interaction time of fast and slow 
wind will also change. Since we suspect that stellar wind 
influences the expansion velocities of the LBV nebulae, all 
effects that determine the wind velocity will automatically 
also influence our expansion velocities. In this context we 
will only briefly discuss one process that changes the stel- 
lar wind velocity, namely metallicity or more general abun- 
dances, since e.g. rotation might change the abundance in 
the outer layers of the stars. If we have stars with lower 
metallicity/abundances — such as in the LMC objects — they 
generally have slower wind velocities, at least in their main- 
sequence phase. Nothing is known so far about the wind 
velocity of stars in the LBV phase with lower metallic- 
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Table 2. Comparison of the parameters of Galactic and LMC LBV nebulae. Numbers separated by slashes indicate that the 
nebula consists of two parts. The sizes are given as diameters or the minimum and maximum extension (smallest and largest 
axes). 
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ity/abundances. Therefore, different abundances — initially 
or, e.g., due to mixing from rotation — will also influence 
the expansion velocities. 

- Eruptions: For some LBVs (e.g. r\ Car, P Cyg) we know 
that they underwent a giant eruption, and at least parts of 
the nebulae were created in these outbursts. Whether or not 
all LBVs had such an outburst is not known (due to the 
lack of historical records). Therefore, the expansion veloc- 
ities of the nebulae might be different since some resulted 
from outburst while outer nebulae formed through interac- 
tion. If we assume that all LBV nebulae were created in 
outbursts, then their outbursts can easily have a different 
strength and duration, which would also result in a variety 
of expansion velocities. Since we do not know what the un- 
derlying mechanism of the giant eruption is, this process 
could be highly dependent on metallicity and then explain 
why LMC LBVs expand so much more slowly. 

- Environment: Whatever created the nebulae, wind or erup- 
tion, the surrounding ISM will also have an impact on the 
expanding nebulae. Density and temperature of the envi- 
ronment that harbors the LBVs will have a noticeable influ- 
ence on the nebulae. By sweeping up ISM, the nebulae will 
slow down and start mixing. Instabilities will occur (e.g., 
Rayleigh-Taylor, Kelvin-Helmholtz) that can lead as far as 
to the final disruption of the nebula and possible outflows as 
seen in, e.g., S 1 19. So far, the closer environment in which 
the LBVs are situated has barely been studied, so that noth- 
ing is known about its influence on the LBV nebulae. 

Some of the above mentioned parameters — perhaps all — are 
responsible for the different expansion velocities of the LBVs 
in the LMC and the Milky Way. As a natural consequence 
a comparison of the nebula with stellar parameters would be 



made. For a specific reason this comparison has explicitly been 
awoided. LBVs do change their stellar parameters significantly 
in timescales as short as some years. Therefore, a simple com- 
parison of stellar parameters of LBVs with each other and with 
the nebula parameters would lead to invalid conclusions. The 
stellar parameters, especially the wind velocity, were different 
at the time the nebulae formed. Comparing the stellar wind, as 
detected in the current spectra, would not help to pin down the 
different properties of the nebulae in the LMC and the Milky 
Way. 

LBV nebulae in the LMC are similar in shape and size to 
those in the Milky Way, but their expansion velocities are dif- 
ferent. Nebulae are either formed in giant eruptions or due to 
wind-wind interaction. We conclude that the formation of LBV 
nebulae seem similar for each nebula, but must occur with dif- 
ferent strengths to account for the different velocities. At least 
to some degree, metallicity seems to play a role in the forma- 
tion of LBV nebulae as the different expansion velocities of 
LMC and Galactic nebulae indicate. 
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